Different wavelength oscillations in the conductance of 5d metal atomic chains 
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Combining ab-initio and self-consistent parametrized tight-binding calculations we analyze the 
conductance properties of atomic chains of 5d elements like Au, Pt and Ir. It is shown that, in 
addition to the even-odd parity oscillations characteristic of Au, conduction channels associated 
with the almost full d bands in Pt and Ir give rise to longer periods which could be observed in 
sufficiently long chains. The results for short chains are in good agreement with recent experimental 
measurements. 
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A traditionally idealized textbook example, the atomic 
linear chain, has recently become an actual system that 
can be explored experimentally. The formation of atomic 
chains several (up to 7-8) atoms long has been achieved 
in recent years using experimental techniques like scan- 
ning tunneling microscope and mechanically controllable 
break junctions In both cases, evidence has been 
found of chain formation in the last stages of pulling an 
atomic contact for certain metallic elements like Au, Pt 
and Ir. It has been suggested that the stability of the 
chains for these 5d elements as opposed to other met- 
als is related to relativistic effects involving transfer of 
charge between sp and d bands 2] . 

These atomic chains exhibit, on the other hand, in- 
teresting conductance properties. Recent experimental 
results indicate the presence of conductance oscillations 
with period ~ 2a (a being the interatomic distance) after 
averaging over many realizations of the chain Q . While 
in the case of Au the oscillations are superimposed to an 
almost constant background of the order of the quantum 
of conductance Gq = 2e^//i, in the case of Pt the mean 
value of the conductance exhibits a continuous decrease 
from ^ 2.5 to ~ I.OGq. On the other hand, for Ir chains 
the conductance varies between ~ 2.2 and ~ I.SGq with 
a less clear oscillatory behavior. 

Theoretical studies of the conductance of atomic chains 
have been so far mainly restricted to idealized models 
or to ab-initio calculations for light metallic elements like 
Na, Si or C Sj , but realistic calculations for the actual 5d 
metal chains, in which the role of d orbitals is expected 
to be relevant, are still lacking From the theoretical 
point of view even-odd conductance oscillations are al- 
ready present at the level of simple tight-binding (TB) 
models for the half- filled case This is an interfer- 
ence phenomena arising from the commensurability of 
the Fermi wavelength and the lattice spacing. This be- 
havior has also been observed in more realistic calcula- 
tions based on ab-initio methods for Na and C This 
simple explanation can account qualitatively for the be- 
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FIG. 1: Model geometry considered for the calculation of the 
conductance. The electrodes are represented by semi-infinite 
fee crystals grown in the (111) direction. 



havior in the case of Au, characterized by a full 5d band 
and a nearly half-filled 6s band. 

For single atom Pt contacts the value around 2.5 Go 
has been attributed to the contribution of additional 
channels associated with the 5d orbitals j^. Thus the 
observation of oscillations with period 2a is not at all ob- 
vious in this case and deserves a detailed analysis. The 
aim of this work is to understand the possible mecha- 
nisms leading to this behavior. We will show that the 
channels associated with the almost full d bands in Pt 
and Ir exhibit transmission oscillations with longer peri- 
ods which can be associated with the Fermi wave vector 
of the bands in the infinite chain. 

We shall consider model geometries like the one de- 
picted in Fig. n in which the atomic chain is connected 
to bulk electrodes represented by two semi-infinite fee 
perfect crystals grown along the (111) direction. The 
interatomic distance a will be allowed to take larger val- 
ues than in the relaxed structure of the infinite chain to 
account for the effect of the applied stress. For simplic- 
ity we shall concentrate first in the analysis of the ideal 
linear geometry discussing in a second step the effect of 
possible structural deformations. 

A convenient approach to quantum transport calcu- 
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FIG. 2: Band structure of the infinite Pt chain for two dif- 
ferent interatomic distances. The symbols and the full lines 
correspond to the ab-initio and to the fitted TB bands respec- 
tively. The bands are classified by the quantum number m 
corresponding to the projection of the angular momentum on 
the chain axis. The arrows indicate the crossing of the Fermi 
level for the m = and the m — ±1 bands. 



lations in atomic size conductors uses a self-consistent 
parametrized TB model in which the main features of 
the bulk band structure are included 0- This type of 
calculation has been fruitful to obtain information about 
the conduction channels in one atom contacts . In the 
particular case of atomic chains a natural starting point 
is the band structure of the infinite linear chain whose 
main properties are discussed below. 

The electronic structure and total energy calculations 
for the infinite chains were performed using the WIEN97 
code This is an implementation of the linearized 

augmented plane wave method based on density func- 
tional theory. For exchange and correlation we have used 
the local density approximation (LDA), with the corre- 
lation part as given by Perdew and Wang [l^ . 

We will discuss the band structure of the chains taking 
the case of Pt as a representative example (Au and Ir ex- 
hibit a rather similar band structure although displaced 
with respect to Ep)- Fig. |21shows the bands around the 
Fermi energy for Pt at two different situations with in- 
creasing stress: one close to the calculated equilibrium 
situation (a = 2.33 A) and the other corresponding ap- 
proximately to the maximum stress before breaking the 



chain (a = 2.90 A). The fits of these bands using a Slater- 
Koster TB scheme are shown as full lines. Several fea- 
tures of this band structure are worth commenting: 1) 
Symmetry considerations allow to classify the bands ac- 
cording to the projection of the angular momentum along 
the chain axis m (we take the z axis along the chain di- 
rection). 2) Close to the Fermi level there is an almost 
flat filled two- fold degenerate band with d^y and dx2_y2 
[m = ±2) character. The other partially filled and more 
dispersing bands have s —pz — d^^ (m = 0) and Px — d^z 
or py — dyz [m = ±1) character (see labels in Fig. [21). 3) 
When the chain is elongated the dispersion of the bands is 
reduced and a net charge transfer from s to d is observed, 
tending rather fast to the "atomic" charge configuration 
5(i^6s^. As a consequence the upper m — band with a 
predominant s character gets closer to the half-filled sit- 
uation [kpa ^ 7i'/2) which, as commented above, would 
give rise to even-odd conductance oscillations. Notice on 
the other hand that the m = ±1 bands cross the Fermi 
level at a lower wave vector. 

In order to describe the electronic properties for a 
model geometry like the one of Fig. ^ the TB Hamil- 
tonian is written a,s H = J2 f^hj,a,pcl^a,a'^j,l^,cr = H^ham 
+ Vl + Vr + Hl + Hr, where Hchazn and Hl,r describe 
the uncoupled chain and the left and right electrodes re- 
spectively; Vl.r being a coupling term between the chain 
and the electrodes which, in the model geometry of Fig. 
n corresponds to the hopping elements connecting the 
outermost atoms of the chain and their three nearest 
neighbors on the electrodes surface. The matrix elements 
hi,j,a,f3, where z,j and a, (3 design sites and orbitals re- 
spectively, are taken from the fits to the chain and to 
the bulk ab-initio band structure. As a self-consistency 
condition we impose local charge neutrality by fixing the 
charge of the d orbitals and sp orbitals within the chain 
equal to the ones of the infinite chain The cor- 

responding condition is also imposed on the electrodes 
surfaces. 

Once the TB Hamiltonian has been built the conduc- 
tance G{E) = GqTt [i''{E)i{E)] of the chain is calculated 
in terms of the matrix elements of the Green function op- 



erator G'^(E) ~ lim^ 



E- 



H 



using 



i{E) = 2tI^'^{E)GI 



N 



(1) 



where TL,Fi. are the matrix tunneling rates connecting 
the chain to the leads |^. The transmission matrix iH 
can then be diagonalized in order to obtain the transmis- 
sion eigenchannels and eigenvalues t„ which characterize 
completely the linear transport properties of the system. 

In order to illustrate the main features and the dif- 
ferences between Au, Pt and Ir, we show in Fig. 13 the 
density of states and the energy dependent transmissions 
for a iV = 5 chain of these metals at an intermediate 
elongation. As has been shown in previous works |^ , Au 
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FIG. 3: Local density of states (LDOS) at the central atom 
and total transmission for An, Pt and Ir chains with N — 5 
at an intermediate elongation. The LDOS is decomposed in 
s (full line) d (dotted line) and p dashed line with the same 
normalization in the three cases. 



chains are characterized by a single conducting channel 
around the Fermi energy with predominant s character. 
The transmission of this channel lies close to one and 
exhibits small oscillations as a function of energy resem- 
bling the behavior in a single band TB model 0. For 
odd N charge neutrality provides a strong mechanism 
for the almost perfect quantization of the conductance, 
as discussed in Ref. Q for one atom. 

In the case of Pt the contribution from the almost filled 
5d bands becomes important for the electronic properties 
at the Fermi energy. There are three conduction chan- 
nels with significant transmission at Ep: one due to the 
hibridization of s — and d^^ orbitals, and other two 
almost degenerate with — dxz and py — dyz character 
respectively. One can naturally associate this channel 
distribution with the band structure presented in Fig. |5J 

The contribution of the 5d orbitals is even more impor- 
tant in the case of Ir where a fourth channel related to 
the lower m = band exhibits a significant transmission. 

Let us analyze the evolution of the conductance with 
the number of atoms in the chain. We find that this 
evolution is rather sensitive to the elongation, specially 
in the case of Pt and Ir (for Au the conductance exhibits 
small amplitude even-odd oscillations, ~ 0.04Go, which 
remain basically unaffected upon stretching). The center 




FIG. 4: Evolution of the conductance with A'^ for different 
values of the interatomic distance a. 



panel in Fig. 21 shows the total conductance of Pt chains 
as a function of N for different interatomic distances. 
Notice that a slight increase of a from the equilibrium 
value leads to a global decrease of the conductance with 
N together with the appearance of even-odd oscillations. 
For a = 2.9 A the total conductance drops from 2.46*0 to 
around IGq for N > A, while the oscillations amplitude 
is of the order of 0.2Go. 

A deeper understanding of this behavior can be gath- 
ered when analyzing the evolution of the conductance 
and its channel decomposition for even longer chains. 
This is illustrated in Fig. [51 for the case of Pt and Ir. As 
can be observed the decrease of the total conductance of 
Pt for iV < 7 — 8 corresponds actually to a long period 
oscillation in the transmission of the two nearly degen- 
erate channels associated with the to = ±1 bands. This 
period can be related to the small Fermi wave vector of 
these almost filled d bands, as indicated by the left arrows 
in Fig.O In addition, the upper to = band crossing the 
Fermi level is close to half-filling giving rise to the even- 
odd oscillatory behavior observed in the transmission of 
the channel with predominant s character. The lower 
m — band tends to be completely filled as the chain is 
elongated and the corresponding channel is nearly closed 
for short chains. However, for intermediate elongations 
one can still appreciate a very long period oscillation in 
its transmission, rising up to ~ 0.5Go, as can be observed 
in the upper panel of Fig. [SI 
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FIG. 5: Channel decomposition of the conductance for Pt 
and Ir as a function A'^ for intermediate elongations. The 
legends indicate the symmetry of the associated bands in the 
infinite chain. The dotted lines in the case of Ir correspond to 
functions of the form ~ coa^ {kF,iNa) where kF,i is the Fermi 
wave vector of the associated bands, which help to visualize 
the oscillatory behavior. 



In the case of Ir both the m = ±1 and the lower m = 
are less filled than in Pt leading to four partially open 
channels. The corresponding values of the Fermi wave 
vector fix its overall oscillatory behavior as illustrated in 
the lower panel of Fig. [S] 

The results presented so far correspond to the idealized 
geometry depicted in Fig. ^ Deviations from this situa- 
tion can arise both from disorder in the atomic positions 
and possible structural distortions. Previous studies indi- 
cate that the relaxed structure for Au chains corresponds 
to a zig-zag configuration [isj . Our ab-initio calcula- 
tions indicate that this is also the case for Pt and Ir, the 
main effect of the external stress being first to suppress 
the zig-zag deformation and then to increase the inter- 
atomic distances. This deformation has some effect on 
the band structure, mainly consisting in the breaking of 
the degeneracies present in the linear case. Although this 
would affect the results close to the relaxed situation, we 
have checked that the main conclusions regarding the ap- 
pearance of long period oscillations associated with the 
d bands are still valid. On the other hand, molecular dy- 
namics simulations suggest a non-uniform distribution of 
bond lengths along the chain We have found that 

this type of deformation as well as the inclusion of mod- 



erate disorder in the atomic positions do not give rise to 
qualitative changes in the conductance behavior as long 
as they remain within the ranges predicted in Ref . |Q| . 

In conclusion we have analyzed the conductance prop- 
erties of Au, Pt and Ir atomic chains. In addition to the 
well understood even-odd parity effects of Au, in metals 
like Pt and Ir the partially filled 5d bands give rise to 
a more complex oscillatory pattern in which more than 
a single wavelength can be identified. As a general rule 
we predict that the transmission corresponding to each 
conduction channel oscillates as ~ cos'^{kF,iNa) where 
kp^i is the Fermi wave vector of the associated band in 
the infinite chain. In the case of Pt the charge transfer 
from s to d orbitals that takes place upon stretching leads 
to a nearly half-filled s band with the consequent even- 
odd oscillations in the corresponding channel, while the 
channels associated with the almost filled d bands ex- 
hibit longer period oscillations. The total conductance 
for short chains {N < 7 — 8) exhibits an overall de- 
crease with superimposed even-odd oscillations in qual- 
itative agreement with the experimental results of Ref. 
. In the case of Ir we find four partially open channels 
with rather similar wavelengths leading to a less clear 
oscillatory pattern. Our predictions for the channel de- 
composition could be checked experimentally by similar 
techniques as the ones used in Refs. [lO, [Tgl • Advances 
in the fabrication techniques would allow to test these 
predictions for longer chains. 
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